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The uncertainty principle generally prohibits determination of certain pairs of 
quantum mechanical observables with arbitrary precision and forms the basis of in- 
determinacy in quantum mechanics [ij, Ql- It was Heisenberg who used the famous 
gamma-ray microscope thought experiment to illustrate this indeterminacy [3]. A 
lower bound was set for the product of the measurement error of an observable 
and the disturbance caused by the measurement. Later on, the uncertainty relation 
was reformulated in terms of standard deviations [4, 5], which focuses solely on in- 
determinacy of predictions and neglects unavoidable recoil in measuring devices [6|. 
A correct formulation of the error-disturbance relation, taking recoil into account, 
is essential for a deeper understanding of the uncertainty principle. However, the 
validity of Heisenberg's original error-disturbance uncertainty relation is justified 
only under limited circumstances ppJo']. Another error-disturbance relation, derived 
by rigorous and general theoretical treatments of quantum measurements, is sup- 
posed to be universally valid[10l 111] . Here, we report a neutron optical experiment 
that records the error of a spin-component measurement as well as the disturbance 
caused on another spin-component measurement. The results confirm that both 
error and disturbance completely obey the new, more general relation but violate 
the old one in a wide range of an experimental parameter. 

The uncertainty relation was first proposed by Heisenberg[3| in 1927 as a limitation of simul- 
taneous measurements of canonically conjugate variables due to the back action of the measure- 
ment: the measurement of the position Q of the electron with the error e{Q), or "the mean error", 
induces the disturbance ri{P), or "the discontinuous change", of the momentum P so that they 
always satisfy the relation 

<Q)v{P) - \. (1) 

where Ti is Planck's constant divided by 27r (here, we use ^ for consistency with modern treat- 
ments). In a mathematical derivation of the above relation from the commutation relation 
QP — PQ = ih, Heisenberg [Sj used the reciprocal relation a{Q)a{P) > | for standard deviations 
cr{Q), (t{P) of position and momentum, which was proved shortly afterwards by KennardQ for 
arbitrary wave functions. This relation was generalized to arbitrary pairs of observables A, B hy 
Robertson j5(| as 

a{A)a{B)>\m[A,B]m (2) 

in any states ji/') with a{A),a{B) < oo. Here, [A, B] represents the commutator = 
AB — BA and the standard deviation is defined as cr(A)^ = {■)p\A'^\ip) — {tp\A\'>p)'^ . Robertson's 
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relation (Eq. [5]) for standard deviations has been confirmed by many different experiments. In 
a single-sHt diffraction experiment [l2j the uncertainty relation, as expressed in Eq. [51 has been 
confirmed. The slit width determines the position spread while the diffraction pattern on the 
screen shows the momentum distribution: the pattern for a narrower slit gets wider and vise versa. 
A trade-off relation appears in squeezing coherent states of radiation fields Starting with 

a theoretical proposal [l5!| and the first experimental generation of squeezed states [l^, many 
experimental demonstrations have been carried out (TtI. Il8l| . 

Robertson's relation (Eq. ^ has a mathematical basis, but has no immediate implications 
to limitations on measurements. This relation is naturally understood as limitations on state 
preparation or limitations on prediction from the past. On the other hand, the proof of the 
reciprocal relation for the error e(^) of an ^-measurement and the disturbance ri{B) on observable 
B caused by the measurement, in a general form of Heisenberg's error-disturbance relation 



is not straightforward, since Heisenberg's proof used an unsupported assumption on the state just 
after the measurement (lOj . Recently, rigorous and general theoretical-treatments of quantum 
measurements have revealed the failure of Heisenberg's relation (Eq. [T|), and derived a new 



Here, the error e{A) is defined as the root-mean-square (rms) deviation of the output operator Oa 
actually measured from the observable A to be measured, whereas the disturbance r](B) is defined 
as the root-mean-square of the change in the observable B during the measurement fldi (see 
Methods A for details). The additional second and third terms result mathematically from non- 
commutability between B (A) and the error (disturbance) operator (Eq. 235 in Ref. [11] ) . In par- 
ticular, they imply a new accuracy limitation e(^) > ^ |(?/;| _B] |'(/') | <t{B)^^ for non-disturbing 
{'i]{B) — 0) measurements and a new disturbance limitation ry(i?) > ^ |('0|[^, cr{A)^^ for 

noise-free (e(j4) = 0) measurements, instead of e{A) ~ cx) or r/(i?) ~ oo as derived from the 
Heisenberg type relation (Eq. 2). 

In this letter, the universally valid error-disturbance relation (Eq. |4|) is experimentally tested 
for neutron spin- measurements [2^. We set A and B as the x- and y-component of the neu- 
tron i-spin. (For simplicity, ^ is omitted for observables of each spin component.) The error 
e{A) and the disturbance ri{B) are defined for a measuring apparatus called Ml, so that the 
apparatus Ml measures the observable A = ax with error e{A) and disturbs the observable 
B ~ ay with disturbance 'r]{B) during the measurement. To control the error e{A) and the dis- 
turbance Tj{B), the apparatus Ml is designed to actually carry out the projective measurement 
of Oa = ct^ = CTj; cos (f) + ay sin instead of exactly measuring A = a^ hy detuning the azimuthal 
angle of ct^, which is an experimentally controlled parameter, so that e{A) and ri{B) of are 
determined as a function of (j) (see Methods B for details). Since the output operator Oa and 
the observable A to be measured are not simultaneously measurable, their difference is not a 
directly detectable quantity, and likewise the same is true for the change in the observable B 
during the measurement. On this ground, the notions of the error e{A) and the disturbance 
ri{B) have been often claimed to be experimentally inaccessible [2ll [22j . In order to overcome the 
alleged experimental inaccessibility, we follow the theoretical analysis (p. 387 in Ref. [llj]) which 
proposes a method to determine the error t{A) and the disturbance rj{B) from experimentally 
available data. A different proposal for the experimental demonstration of the same relation 
(Eq. 3]), which exploits the weak-measurement technique, was recently published by Lund and 
Wiseman |23i]. 



e[A)^{B)>-\mA.Bmi 



(3) 




e{AUB)+e{A)a{B)+a{A)^{B) > -MiA^BMl 



(4) 
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The error e{A) is deterrained by the data from the apparatus Ml, and the disturbance ri{B) 
is determined by the data from another apparatus called M2 which carries out the projective 
measurement of B on the state just after the Ml-measurement. Thus, the experiment is based 
on the successive projective measurements of two non-commuting observables Oa in Ml and B 
in M2 as depicted in Fig. 1. 
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FIG. 1: Experimental concept for the demonstration of the error-disturbance uncertainty relation in 
successive measurements: After preparing an initial state \tp), apparatus Ml is assumed to measure an 
observable A (red region). The error e{A) in the ^-measurement is experimentally controlled by detuning 
apparatus Ml to measure Oa instead of A (light red). After this measurement, the state is projected 
onto one of the eigenstates of Oa, inevitabely influencing the subsequent measurement of observable B, 
performed by apparatus M2 (yellow region). The disturbance ry(-B) on the _B-measurement depends on 
the detuning of Ml, i.e. the error of the yl-measurement. The successive spin-i measurements of Oa 

and B result in four possible outcomes, denoted as (-f +), (H — ), ( — h) and ( ), from which error 

e{A) and disturbance ri{B) are quantitatively determined. 



For these measurements, the neutron beam first passes the preparation stage of the initial 
state \tp). The apparatus Ml has two possible outcomes, i.e. -|-1 and —1, corresponding to 
measurement operators E'^{±1) = i(J±cr^). The disturbance on the observable B caused by the 
apparatus Ml is detected by the apparatus M2, which also yields either -1-1 or —1, corresponding 
to measurement operators i?^(±l) ^(^i ay). Thus, the successive measurements carried out 

by Ml and M2 finally result in four intensities denoted as (++), (H — ), ( — h), and ( ). The 

setup of the neutron spin experiment is depicted in Fig. 2. The azimuthal angle (/) of tr^ is tuned 
between and -I in a way that a trade-off for the error and the disturbance occurs. As seen from 
Ea.(|lip and Eq. (jl3l) error e{A) and disturbance ri{B) are obtained by performing the successive 
measurement of Oa and B on different states. For measuring error and disturbance in the state 
IV') = \+z), the auxiliary states \ — z), \+x) and are prepared likewise (see Methods C for 
more details). 

The experiment was carried out at the research reactor facility TRIGA Mark II of the Vienna 
University of Technology (TU Vienna). The monochromatic neutron beam with a mean wave- 
length of 1.96 A propagates in -|-y-direction. The beam is approximately 99% polarized crossing 
a bent Co-Ti super-mirror array (polarizer) |2^ . Two analyzing super- mirrors (analyzers) are 
adjusted to higher incident angles so that the second order harmonics in the incident beam are 
suppressed. The final intensity was about 90 neutrons/s at a beam cross section of 10 (vertical) 
X 5 (horizontal) mm^. A '^He monitor detector is used for normalization in order to correct 
statistical fluctuations evoked by the reactor power. A BF3 detector with high efficiency (more 
than 99%) is used for the experiment. To avoid unwanted depolarization a static guide field 
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PREPARATION 

FIG. 2: Illustration of the experimental setup for demonstration of the universally valid uncertainty 
relation for error and disturbance in neutron spin-measurements. The neutron optical setup consists of 
three stages: preparation (blue region), apparatus Ml performing measurement of observable Oa = cr(^> 
(red region) and apparatus M2 carrying out measurement of observable B = ay (yellow region). In 
the preparation stage a monochromatic neutron beam is highly (99%) polarized in -l-z-direction by 
passing through a super-mirror spin polarizer. The first DC coil (DC-1) produces a magnetic field in x- 
direction (Bx) which can be used to rotate the inital polarization vector about the x-axis. By additionally 
exploiting Larmor precession about the -f z axis, induced by the static guide field present throughout the 
entire setup, and varying the position of DC-1 arbitrary initial spin states can be produced at end of the 
preparation stage (up to an irrelevant phase factor). The projective measurement of observable Oa is 
realized by similar components: the prepared state rotates about the z-axis due to Larmor precession. 
Hence, by properly placed DC-2 coil the spin component to be measured can be projected towards the 
+z direction, where it is reflected by a super-mirror analyzer (Analyzer-1). After passing through this 
first analyzer in the +z state, DC-3 produces the eigenstate \±(j>) of a^. In the same manner, apparatus 
M2 performs the measurement of observable B on the eigenstate |±(/!>), which constitutes the source of 
the disturbance on the second measurement. The combination of the projective measurements of 
and ay gives four count rates at the neutron detector in the downstream of the beam. The error e{A) 
and the disturbance ri{B), as well as the standard deviations of each measurement a{A) and a{B), are 
determined from the expectation values of the successive measurement. 



pointing in -|-z-direction with a strength of about 10 Gauss permeates rectangular Helmholtz 
coils. In addition, the guide field induces Larmor precession, which, together with four appropri- 
ately placed DC spin rotator coils, allows state preparation and projective measurements of Oa 
in Ml and B in M2. (see Methods D for more details). 

To test the universally valid uncertainty relation stated in Eq. 21 the standard deviations 
(y{A), (t{B), the error e{A) and the disturbance r]{B) are determined. The measurement of the 
standard deviations a-{A) and <t{B) is carried out by Ml and M2 separately, whereas error e{A) 
and disturbance ri{B) are determined by successive projective measurements utilizing Ml and 
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Combined Outcomes of Ml and M2 



FIG. 3: Normalized intensity of the successive measurements carried out by apparatus Ml and M2. The 

successive measurement of Ml and M2 has four outcomes, denoted as (++), (H — ), ( — h) and ( ). 

Intensities, according to the corresponding outcomes, are depicted for each initial spin states, i.e. \+z}, 
\—z), \+x) and \+y)- Three sets for detuning parameter = 0,40 and 90deg are plotted. The error 
e{A) and the disturbance ri{B) are determined from these 16 intensities, for each setting of the detuning 
parameter 0. 



M2. Typical experimental data sets, for miscellaneous detuning angles 0, are depicted in Fig.|3l 
The resulting values of e{A) and r]{B), together with the theoretical predictions e{A) = 2 sin ^ 

and ri{B) — \/2cos(f> (see Methods B), are plotted as a function of the detuning parameter cj) 
in Fig.m The trade-off relation of e{A) and ri{B) is in good agreement with theory: when one 
observable is measured more precisely, the other is more disturbed. The vertical and horizon- 
tal error bars in Fig.|4] contain the statistical fluctuations of the measurement as well as the 
systematical misalignments of coil position and current values resulting in an angle deviation 
of 1.6 degrees. The final results were obtained by taking the contrast (~ 96%) of the entire 
measurement into account. 

From the terms obtained above (error e{A), disturbance ri{B), standard deviations cr{A) and 
cr{B)), the Heisenberg error-disturbance product €{A)r]{B) and the left side of the new relation 
(Eq. SI) e{A)i]{B) + €{A)a{B) + a{A)r]{B) are plotted as a function of the detuned azimuthal 



6 



1.5 




Detuning Angle (j) [rad] 

FIG. 4: The trade-off relation between the error e{A) (blue), and the disturbance r]{B) (green), obtained 
as a function of detuning the azimuthal angle (j). The theory predicts the dependence of (j): e{A) = 2 sin ^ 
and 77(B) = cos (j)- The values at (j) — are €{A) = and 7?(-B) = \/2 whereas the values at 4> — ^ are 
e{A) = y/2 and r]{B) = 0. It is impossible to accurately measure both observables A and B at the same 
time. Vertical and horizontal error bars include statistical and systematical errors. 

angle (j) in the upper panel of Fig.[S] This figure illustrates the fact that the Heisenberg product 
is always below the calculated limit, and that the new sum is always larger than the limit in 
the scanned range of </>. This clearly demonstrates that the Heisenberg-type error-disturbance 
relation (Eq.[3]) with the single product of the error e{A) and disturbance ri{B) is violated whereas 
the new relation (Eq. 01) consisting of three terms is always satisfied. 

The technique utilizing various different incident states for the estimation of the effects of the 
quantum operation reminds us of quantum process tomography [2 5]. Although a parameter for 
the error is experimentally controlled in the experiments, it is easily extendable for uncontrolled 
and fluctuating parameters. Here, we concentrate on the situation where the full trade-off relation 
between error e{A) and disturbance r](B) occurs. It is worth noting that the mean value of the 
observable A is correctly reproduced for any detuning angle (j), i.e., {+z \Oa\ + z) = {+z \A\ + z), 
so that the projective measurement of Oa reproduces the correct probability distribution of A, 
while we can detect the non-zero rms error e{A) for (p 0. What has been accepted by the 
uncertainty principle is the existence of an unavoidable trade-off between measurement accuracy 
and disturbance, but this principle has eluded a satisfactory quantitative description for a long 
time. Our result is the first evidence for the validity of the new relation (Eq. S]) proposed as 
a universally valid error-disturbance relation by one of the present authors^lO, 11], whereas the 
failure of the old relation (Eq. [S]) is also illustrated. Our results witness that the new relation 
solves a long-standing problem of describing the relation between measurement accuracy and 
disturbance. Our demonstration sheds a light not only on fundamental limitations of quantum 
measurements but also on technical limitations of precise measurements such as gravitational- 
wave-detection (26l - [29l | and quantum information processing. 
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FIG. 5: Experimentally determined values of various products consisting of error e, disturbance rj, 
and standard deviations a together with theoretical predictions. (Upper panel) The left side of Eq.Q 
e{A)i]{B) + e{A)a(B) + a{A)i]{B) = 2^sin | cos + 2 sin | + \/2cos(j} (orange) and the Heisenberg 
product e{A)ri{B) = 2^2 sin ^ cos cj> (red) are plotted as a function of the detuned azimuthal angle ip. The 
Heisenberg product is always smaller than the calculated limit ^ | {^/jI [v4, -B] IV') I = ^1 {+z \ [a^ , o-y] \ + z)\ = 
1 (depicted as the dashed line). In contrast, the new sum is always larger than the limit. This is a clear 
evidence of the violation of the old uncertainty relation, in a sense of Heisenberg's error-disturbance 
relation, and the validity of the new relation consisting of the three-terms. (Lower panel) The two 
additional product-terms a{A)ri{B) (green) and e{A)a{B) (blue) in the new relation are plotted together 
with the theoretical predicted curves: €{A)a{B) = 2 sin I and a{A)rj{B) = ^/2 cos (j>. Error bars of the 
experimental data include statistical and systematical errors. 



Methods 



A) Universally valid uncertainty relation 

Any measuring apparatus M is, in principle, modeled by the unitary operator U describing 
the time evolution of the composite system of the measured object S and the probe system P 
during the measuring interaction and the meter observable AI of P actually measured after the 
measuring interaction [11]. If the initial states of the object and the apparatus is and |^), 
respectively, the root- mean-square (rms) error e{A) of M for measuring an observable A of S 
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and the rms disturbance ri{B) of M caused on an observable _B of S is defined as 

eiA) = \\[U\l®M)U-A®I]\i,)\0\\, (5) 
7?(S) = \\[uHB^I)U-B^I]\^)\0\\. (6) 

Then, it is mathematicahy proved [13, [HI that Eq.H] holds for any unitary operator U for S + P, 
observable M of P, and state vector \tp) of S and |^) of P. Suppose that the apparatus M has 
a family {M,n} of the measurement operators [25l|. This means that the measuring apparatus M 
has possible outcomes m with probability p{m) = \\Mm |P and the state of the object S after 
the measurement is Mm \ip) /WMm IV') II- In this case, the rms error and rms disturbance are give 

e{Af = ^||M„,(m-A)|^)f, (7) 
viBf = ^||[Af„,,B]|^)f. (8) 

m 

If {Mjn} consists of orthogonal projections, the measurement is called a projective measurement. 
In this case, Eq. [7] can be simplified as e{A) = \\{0a — A) |-0) || by the Pythagorean theorem, 
where Oa — J2m "'^^■m is called the output operator. 



B) Error and disturbance in spin measurements: Theoretical determination 

In the experiment, we test the universally valid uncertainty relation (Eq. ^ for observables 
A — Ux and B = ay, while the initial state \tp) is \+z) and the measuring apparatus M — Ml 
is considered to carry out the projective measurement of Oa = cTip — cos (pax + siTicjjay. Thus, 
the apparatus Ml is described by measurement operators £''^(+1) = (1 + a^)/2 and E'^{—1) = 
(1 - a^)/2 with Oa = J2x=±i xE'>'{x). From Eqs. ^ and ©, we have 

e{A) = ||(a^-a,)|V')|| =2sin^, (9) 
r]{B) = \/2||[(7^,(7j,]|^/.)|| = \/2cos0. (10) 



C)Experimental determination of error and disturbance 

In the experiment, we determine e(A) and rj{B) from statistically available data obtaine d by 
successive neutron spin-measurements. According to the former theoretical analysis (Ref. |ll| . 
p. 387), the error e{A) is determined by mean values of Oa in three different states as 

t{Af = (VI^'IV') + {MO\\i^) + (^|0a|^) + {MAOAA\i^) - (v-KA + i)Oa{a + im 

= 2 + (V'IOaIV') + {Ai;\OA\A^) - {{A + /)V|Oa|(A + /)V), (11) 

where we have used the following abbreviations: \A'il}) = A|?/;) and \{A + I)'4>) = [A + I) l-^). 
Since the apparatus Ml carries out the projective measurement of Oa, in order to determine 
e{A) for the basic initial state we need only to measure the intensities from the apparatus 
Ml in the three auxiliary incident states of Ml corresponding to |'0) , ^IV') j + I) IV')- The 
expectation values expressed in Ea. ljllD are calculated from the measured intensities, depicted in 
Fig.E via 



(12) 
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In order to detect the disturbance on B caused by the apparatus Ml, the apparatus M2 carries out 
the projective measurement of B in the state just after the Ml-measurement. The modified out- 
put operators of the apparatus M2 for the initial state of Ml is given by Ob = X^a; E'^{x)BE'^{x) 

and = J2x E'*'{x)B'^E'f'{x), which describes the mean and the second moment of the observ- 
able B for the initial state of Ml. Then, from Eqs. (189) and (227) of Rcf. [H] the disturbance 
rj{B) is also determined by mean values oi Ob in three different states as 

in{Bf = (VIS'IV-) + (V'loi'V) + {MOb\^) + {Bi^\OB\Bi^) - {{B + I)^\OB\{B + m 

= 2 + (^|Ob|^) + {Bij\OB\B^) - {{B + I)^\Ob\{B + I)^), (13) 

where the illegitimate notations l-BV') and \{B + are used as before and the expectation 
values are given by 

i++ + 1+- + 1-+ + 1 — 

Thus r}{B) is determined in the same manner as e{A). By the relations jV') = \+z), A jV') = 
B\ii) = i\-z), {A + I)\'il;) ^ V2\+x), and {B + I)\ip) = V2\+y), where we set |+a;) = 
(l+z) + |— z))/\/2 and \+y) = {\+z) + i \—z))/^/2, the required spin-states in Eqs. [TT] and [13] 
are generated by spinor- rotations in the experiment (see Fig|2]); the normalization factors are 
confirmed experimentally in spin-rotation measurement. 



D) Successive neutron spin-measurements 



The spin state of the neutrons is controlled by four DC coil spin-turners. The required incident 
states are prepared by first DC coil (DC-1). DC-1 is switched off for the generation of |+z) 
incident state and the spin is fiipped by DC-1 for |— z) preparation. While the incident spin 
state \+y) is generated by applying a ^ rotation around the x-axis, the incident spin state \+x) 
is produced by additionally moving the position of DC-1 one quarter of Larmor rotation period 
due to the guide field. The apparatus Ml consists of the combination of the spin-turner coils, DC- 
2 and DC-3, the guide field and the spin analyzer 1. The spin analyzer in our experiment performs 
the projective measurement of the spin's +z-component. Instead of rotating the analyzer, the 
neutron's spin-component in the x-y plane is unitarily rotated towards the analysis-direction by 
DC-2 and the guide field. DC-3 finally generates the eigenstate |±(/)) = i5"*(±l) \^) (up to phase 
factor), so that the apparatus Ml performs the projective measurement of Oa to obtain the 
mean values of Oa in Ea. pH) . The detuning of A to Oa is adjusted by shifting DC-2 and DC-3 
likewise towards or away from the analyzer 1. The apparatus M2 consists of DC-4 and analyzer 
2 and performs the projective measurement of B on the state just after the Ml-measurement to 
evaluate the disturbance on B caused by Ml. The coil DC-4 turns the y-component of the spin 
into -|-z-direction by ^ rotation around the x-axis, with which analyzer 2 performs the second 
projective measurement of B] additional spin-rotation towards the y-direction is omitted here, 
since only (spin-insensitive) intensity measurement is performed afterwards. By the disturbance 
caused by the apparatus Ml, the passage until the M2 measurement is described by the output 
operator Ob for the incident state of Ml, and hence we obtain the mean values oi Ob in Ea.([T5|) 
from the apparatus M2. Note that the measurement performed by apparatus M2 is the error-free 
^-measurement on the state just after the measurement carried out by Ml. 
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